Recent measurement of AG/G at COMPASS 
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Abstract. We present a preliminary measurement of the gluon polarization AG/G in the nu- 
cleon, based on the spin asymmetry of quasi-real photoproduction events for which a pair of large 
transverse momentum hadrons is produced. The data were obtained by the COMPASS experi- 
ment at CERN using a 160 GeV polarized muon beam scattered on a large polarized 6 LiD tar- 
get. The preliminary helicity asymmetry for the selected events is Ay/D = 0.002 ±0.019(stat.) ± 
1—5 , 0.003 (exp.syst.). From this value, a leading order analysis based on the model of PYTHIA leads 

■ to the gluon polarization in the nucleon AG/G(x g — 0.095, )J. 2 — 3 GeV 2 ) = 0.024 ±0.089(stat.) ± 

\ 0.057(syst.). This value is consistent with parameterizations obtained from QCD fits to the g\ data, 

with a first moment AG = f AG(x)dx < 1, at the same scale. 
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■^j- . The decomposition of the nucleon spin in the contributions from its constituents has 

been a central topic of investigation in polarized lepton-nucleon scattering in the last 20 
years. The EMC measurement of the proton spin structure [1] has shown that only 20 to 

i/-) ■ 30 % of the proton spin could be attributed to the total quark spin AE, in contrast to the 

60 % expected in the naive quark-parton model. In inclusive lepton-nucleon scattering 
the contribution of the gluon spin AG to the nucleon spin can only be measured indirectly 
by studying the Q 2 dependence of the polarized spin-structure functions in QCD. Fits at 
Qh! next-to-leading order (NLO) provide evaluations of AG which are of the order of 0.50. 

The precision of these fits is strongly limited by the small Q 2 range covered by the data at 
any value of x, a situation resulting from the lack of a polarized lepton-nucleon collider. 
In addition, the shape of AG(x) at a fixed Q 2 used as reference, has to be provided as 
an input parameterization. It varies considerably between the different analyses and is 
only poorly constrained by the results of the fits. A direct measurement of the gluon 
polarization AG(jc) / G(x) can be obtained from the helicity asymmetry of the photon- 
gluon fusion (PGF, y*g — > qq) cross-section, which constitutes an important part of the 
experimental program of COMPASS 

The COMPASS experiment [2] is located at the M2 beam line of the CERN SPS, 
which provides a 160 GeV jti + beam, with a natural polarization of —76 ± 5%. The tar- 
get consists in an upstream and a downstream cell, longitudinally polarized in opposite 
directions. Typical target polarizations of 50.0 ±2.5% are obtained. The forward spec- 
trometer is divided in two stages allowing the reconstruction of the scattered muon and 
of the produced hadrons in broad momentum and angular ranges. The trigger system 
provides efficient tagging down to Q 2 = 0.002 GeV 2 . 

The present analysis focuses on the data collected during the 2002 and 2003 runs. 
We only consider quasi-real photoproduction events (Q 2 < 1 GeV 2 ), in which at least 
two charged hadrons are associated to the primary vertex in addition to the incident 
and the scattered muons. The fraction of PGF is enhanced by asking the two leading 



hadrons to have a large transverse momentum: p T l > 0.7 GeV, pj 2 > 0.7 GeV and 
Or) 2 + (Pt 2 ) 2 > 2 - 5 GeV 2 . In total, around 350,000 events are selected. For this high 
pr sample, the measured helicity asymmetry (defined as in [3]) is 

An 

= 0.002±0M9(stat)±0.003(exp.syst), (1) 

where the quoted systematic error accounts for the false asymmetries related to the 
apparatus. Other sources of systematic errors, including the error on the beam and 
target polarizations, are only a few percents of the (small) measured asymmetry, and 
are therefore negligible. The 2004 data are currently under analysis, and represent the 
same amount of data as 2002 and 2003 altogether. 

A Monte-Carlo simulation is needed to extract the gluon polarization from the high p T 
asymmetry. The selected sample of high p T events covers the transition region ranging 
from photoproduction (Q 2 w 0) up to DIS (Q 2 ~ 1 GeV 2 ). For this reason, we chose 
PYTHIA as an event generator because it provides a model for the lepton-nucleon 
interactions [4] at low Q 2 . Two different kinds of processes are generated. In the so- 
called direct processes, the virtual photon takes part in the hard partonic interaction. In 
the resolved processes, it fluctuates to a hadronic state, from which a parton is extracted. 
This parton then interacts with a parton from the nucleon. At Q 2 < 1, the resolved 
processes constitute half of the high p T sample. Their contribution falls to about 10% for 
Q 2 > 1 GeV 2 and becomes negligible for Q 2 > 2 GeV 2 . Note that the analysis requires 
the factorization between the hard and soft parts of the reaction, hence the presence of 
a hard scale /i 2 . As Q 2 < 1 GeV 2 , the scale is provided by the transverse momentum 
of the partons involved in the reaction. Events for which no hard scale can be found are 
classified as low pj. 

The generated events are tracked through a GEANT description of the COMPASS 
spectrometer, and processed using the same reconstruction program as for real data. 
Then, the Monte-Carlo sample of high p T events is selected through the same cuts. 

Only one parameter of PYTHIA had to be changed to reach a good agreement with 
the data: the width of the Gaussian distribution of intrinsic transverse momentum of 
partons within the resolved virtual photon (PARP(99)) was decreased from 1 GeV/c to 
0.5 GeV/c. Fig. 1 presents a comparison between the simulated and real data samples 
of high p T events. Fig. 2 shows how the Monte-Carlo sample of high p T events divides 
into the various PYTHIA subprocesses. The high pj asymmetry can be written in terms 
of the contributions of the different processes: 
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Note that we have neglected the small contributions of the leading process y*g — > q and 
of the low pj scattering events because there is no hard scale allowing a perturbative 
treatment of these subprocesses (low transverse momentum, Q 2 < 1 GeV 2 events). As 
can be seen on Fig. 2, the fraction of photon-gluon fusion events is Rpgf = 0.31. The 



analyzing power a£^ F is the helicity asymmetry of the fig — > fi'qq scattering cross- 
section, d^ F = dAcfpQf-/ dOpQp. It is calculated from the kinematic variables of the 
partonic reaction for each PGF event in the high pj Monte-Carlo sample. Averaging over 
the PGF events, we obtain (aff F /D) = -0.933. The contribution of the PGF process 

to the high pj asymmetry is thus —0.292 x ^f. The contribution of the QCD Compton 
events is calculated in the same way to be 0.0063, using a fit on the world data for the 
virtual-photon deuteron asymmetry Af . 

Resolved photon subprocesses involve either a quark or a gluon from the nucleon. In 
the latter case, they are sensitive to the gluon polarization AG/G, and contribute to the 

signal. The analyzing powers d L J L are calculated in pQCD at leading order [5], and are 
positive for all relevant channels. The polarizations (Af/f) d of the u, d and s quarks 
in the deuteron are calculated using the unpolarized parton distribution functions from 
GRV98, and the polarized parton distribution functions from GRSV2000 [6], at leading 
order. The polarizations of quarks and gluons in the virtual photon (Af/fy are unknown 
as the polarized PDFs of the virtual photon have not yet been measured. Nevertheless, 
theoretical considerations provide a minimum and a maximum value for each A/ r [7], 
called the minimal and maximal scenarios. The total contribution of the resolved photon 
processes to the high pr asymmetry, which ranges between 0.000 + 0.012 x AG/G and 
0.002 + 0.078 x AG I G, will be taken into account in the systematic error on AG/ G. 

The tuning of the PYTHIA parameters relevant to this analysis is an important source 
of systematic errors. This error was estimated by scanning these parameters indepen- 
dently over a range where the agreement between the simulation and real data remains 
reasonable. This resulted in several values for AG/G, all based on the same high pj 
asymmetry, Eq. (1). The value of AG/G appears to depend predominantly on the width 
of the intrinsic transverse momentum distribution for the partons in the photon. For in- 
stance, varying this parameter between 0.1 and 1 GeV/c results in a 30% variation of the 
fraction of photon-gluon fusion Rpcf - NLO effects seem to be small: varying the scale 
and (de) activating parton showers does not affect the result. 
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FIGURE 1. The upper part of these plots show a comparison between the simulated (histogram) and 
real data (points) samples of high pj events, normalized to the number of events. The lower part shows 
the corresponding data/simulation ratio, p (pj) is the total (transverse) momentum of the leading hadron. 
A similar agreement is obtained for the next-to-leading hadron. 




FIGURE 2. Contribution of each PYTHIA subprocess to the Monte-Carlo sample of high pj events. 
On the left: direct processes (photon-gluon fusion, QCD Compton, and leading process); on the right: 
resolved processes; note that qq — > q'q' , qq — > gg and gg — » qq are neglected, as they altogether represent 
0.6% of the sample. 



Using Eq. (2) to extract the gluon polarization from the high p T asymmetry, Eq. (1), 
we finally get: 

AG -> 

— (x g = 0.095, \l 2 = 3 GeV 2 ) = 0.024 ± 0M9(stat.) ±0.057 (syst.). (3) 
G 

This result was compared to the recent distributions of AG(x) / G(x) from AAC [8] and 
LSS [9]. These two distributions, which strongly differ in shape, are almost equal at 
x g = 0.095 and are compatible with our result within 1.5 c The first moments of AG are 
equal to 0.8 ±0.56 and 1.1 ±0.52 for the AAC and LSS fits, respectively. The fraction 
of PGF can also be enhanced by selecting events with charmed hadrons (D and £>*), 
instead of high pj hadrons. The results obtained with this method, which is much less 
model-dependent but suffers from low statistics, will be produced soon. 
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